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ABSTRACT 

We report the detection of velocity variations in the pulsating sdB star, PG 
1605+072. Oscillations are detected at the same frequencies found from photom- 
etry and have amplitudes of up to 14kms _1 for H/3. The strongest oscillation 
found in previous photometric observations is not evident in our spectroscopy or 
photometry, and may be absent due to beating of closely spaced modes. Phase 
differences between spectroscopy and B magnitude photometry imply that max- 
imum brightness occurs not long after maximum radius. We have also found 
evidence of variation in the observed amplitudes of five Balmer lines, with a de- 
crease in amplitude of the strongest mode blueward from H/3. This effect is not 
expected and a longer time-series will be needed to clarify it. 

Subject headings: stars: interiors — stars: oscillations — subdwarfs 

1. Introduction 

The discovery of a class of pulsating hot subdwarfs (sdBs; Kilkenny et al. 1997) has 
created an interest in asteroseismology of evolved stars. Subdwarfs are believed to lie on 
the Extreme Horizontal Branch (Saffer et al. 1994), although their formation and evolution 
is still uncertain. Pulsations in sdBs were predicted theoretically by Charpinet et al. 
(1996) just prior to their discovery, and about 20 pulsators have been discovered to date 
using time-series photometry (see Koen et al. 1999 for a review). The stars (also known 

1 Based on observations made with the Danish 1.54 m telescope at ESO, La Silla, Chile, 
from the South African Astronomical Observatory (SAAO), and from Mt Stromlo Observa- 
tory, Australia. 
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as EC14026 stars, after the prototype) are believed to pulsate in p modes and typically 
have periods 100-200 s and semi-amplitudes <10mmag. Models indicate an opacity bump 
associated with iron ionisation may be the driving mechanism of the pulsations (Charpinet 
et al. 1996, 1997), although recent studies by Koen et al. (1999) find no well-defined 
instability strip in the log g — T c s plane. 

Combining time-series photometry with time-series spectroscopy should help to give 
a better understanding of the pulsations in these stars. We have shown the feasibility 
of time-series spectroscopy with observations of PG1605+072 on a 1.5 m telescope over 
2 hrs (O'Toole et al. 2000). Jeffery & Pollacco (2000) made similar observations of 
KPD2109+4401 (5 hrs) and PB8783 (6 hrs) using the 4.2m William Herschel Telescope. In 
all three cases, velocity variations were detected at the published photometric frequencies. 
However, sdB stars are multimode pulsators and longer time series are needed to resolve 
the individual modes and measure their amplitudes. We present here the first extended 
time-series spectroscopy of a pulsating sdB star. 

Our target was PG1605+072 which has the most extreme properties of all the 
pulsating sdBs detected so far, with high amplitudes (up to ~25 mmag) and long periods 
(~500s). Because of its low gravity (log g ~5.25), this star is believed to have evolved off 
the Horizontal Branch and may oscillate in g, as well as p modes (Kilkenny et al. 1999). 
The power spectrum is complex, with up to 55 indentifiable frequencies. PG1605+072 
displays considerable rotation (-usim = 39kms _1 ), possibly leading to unequally spaced 
multiplet components (Heber et al. 1999). Unlike several of the pulsating sdBs, multicolour 
photometry gives no indication of a cool companion (Koen et al. 1998). 
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2. Observations 

We obtained medium-resolution spectra of PG1605+072 using the DFOSC spectrograph 
on the Danish 1.54 m telescope at La Silla, Chile and the coude spectrograph (A grating) 
on the 74-inch (1.88 m) telescope at Mt Stromlo, Australia. The observations were made 
on 7 nights over an 11-day period in July and August, 1999 (see Table 1). To supplement 
the spectroscopy, we obtained time-series photometry at the Sutherland site of the South 
African Astronomical Observatory (see Table 2). 

EDITOR: PLACE TABLE 1 HERE. 

EDITOR: PLACE TABLE 2 HERE. 

The La Silla data consisted of single-order spectra projected onto a 2K LORAL 
CCD - pixel-binning (to reduce R.O.N.) and windowing (to reduce readout time) gave 
66 x 500 pixel spectra with a total wavelength range of 3700-5000 A and a dispersion of 
1.65 Apixel -1 . The resolution was 2.8 A, set by a slit width of 1.5 arcsec. The exposure 
time was 46 s, with a dead time of about 15 s. The average number of photons per A in 
each spectrum was about 2200. 

The Mt Stromlo data consisted of single-order spectra projected onto a SITe 2Kx4K 
CCD, with a wavelength range of 3800-5000 A. The slit width was ~4 arcsec and the 
dispersion was 0.60 A pixels -1 , leading to a resolution similar to that of the La Silla data. 
The exposure time was 50 s with a dead time of around 25 s. 

The photometry from Sutherland was taken using the Modular Photometer with 
a Johnson B filter and exposure times of 20 s. The GaAs photomultiplier tube in this 
photometer has a good red sensitivity, and using the blue filter causes the passband to 
resemble that of a blue-sensitive photomultiplier tube (O'Donoghue et al. 1998). 
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3. Reductions 

Standard methods for bias subtraction, flat-fielding and background scattered-light 
subtraction were used. Non-linearities on the CCD image were only important in the 
Stromlo data, and were corrected immediately after bias subtraction (Baldry et al. 1998, 
see). One-dimensional spectra were extracted using the Optimal Extraction method (Home 
1986), a variance weighting system for CCD columns. A 3rd order polynomial was fitted 
to the continuum level in both data sets, and the spectra were normalised to a continuum 
value of unity. 

The spectrum of PG1605+072 is dominated by Balmer lines - see e.g. Figure 9 of Koen 
et al. (1998). A cross-correlation technique was used to determine the Doppler shift for five 
Balmer lines (H/3, H7, H£, He and H8), relative to a template spectrum (the average of 20 
high-quality spectra). The spectra and template were prepared in the manner described by 
Baldry et al. (1998). It should be noted that the time-series for H7 does not contain any 
Mt Stromlo data, because that line was contaminated by a bad CCD column. 

Several jumps were found in each time-series, particularly in the La Silla data, 
occurring every 30-40 observations, corresponding to breaks where the He+Ne calibration 
spectra were taken. We decided to process each of these groups separately with a template 
spectrum created for each group. In removing the jumps, we have effectively high-pass 
filtered the dataset, also removing any instrumental drifts present. 

EDITOR: PLACE FIGURE 1 HERE. 

El 



NOTE TO EDITOR: rawvelocity.ps should be over two columns 
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The photometric data were reduced by subtracting the sky background, correcting 
for atmospheric extinction and normalising to the mean intensity of the run, following a 
method similar to O'Donoghue et al. (1997, 1998). Slow drifts were not corrected for in the 
photometry. 

The H/9 velocity curves and B magnitude light curves for each night of observations are 
shown in Figure 1. Oscillations are not as evident in our velocity curve as they are in the 
light curve, since the signal-to-noise per measurement is lower. 

As is common with time-series measurements, the quality of the data varies considerably 
through the data set. We therefore performed frequency analysis using a weighted Fourier 
Transform, where weights were assigned according to the local rms scatter (Kjeldsen & 
Frandsen 1992). 



4. Results and Discussion 

The amplitude spectra of H/3 and H7 are shown in the top two panels on the left of 
Figure 2. The white noise level is 1.07 km s^ 1 for H/3 and 0.93 km s -1 for H7. Oscillations 
are clearly visible at 2102 //Hz and 2743 //Hz (2742.72 //Hz in Kilkenny et al. 1999) in both 
lines but at 1891 //Hz (1891.42 //Hz) only in H7. It is unclear whether this last effect is due 
to the low amplitude of the mode in H/3 due to beating between closely spaced modes. The 
bottom two panels show our B magnitude and white light amplitude spectra (from a 15-day 
multisite campaign by Kilkenny et al. (1999) in May, 1997; white light observations were 
obtained using a CuS0 4 filter) respectively, and are included for comparison. The three 
modes previously mentioned are all present in the B magnitude spectrum. We believe that 
the highest amplitude peak (14kms _1 for H/?) at 2102 //Hz may be a combination of the 
2101.65 //Hz and 2103.28 //Hz frequencies found by Kilkenny et al. - a longer time-series is 
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needed to clarify this. 

EDITOR: PLACE FIGURE 2 HERE. 

E2 

The panels on the right of Figure 2 show the 2050— 2150 juHz region in greater detail. 
The highest amplitude mode found by Kilkenny et al. at 2075.76 /xHz is not present in any 
of our observations. It is not clear whether this is due to variation in amplitudes over the 
2.2 years between observations, or to beating between very closely spaced modes. 

The B magnitude and white light amplitudes differ by around 40%, which is to be 
expected in such hot, blue stars, considering that the effective wavelength of white light is 
redder than the B band. If oscillation phase is wavelength dependent, then this would also 
reduce the amplitude in white light (which is a much broader "band" than B). 

4.1. Oscillation Amplitudes 

In this Letter we present a qualitative comparison of luminosity and velocity 
amplitudes. All amplitudes were measured by the height of their peak. A more detailed 
quantative study will be the subject of a subsequent paper. 

Kjeldsen & Bedding (1995) derived a relationship between luminosity amplitude at 
a given wavelength and velocity amplitude for classical and solar-like oscillations. Using 
their scaling law, which assumes a black body spectrum and adiabatic oscillations, we 
find that the velocity amplitudes in PG1605+072 are 2-3 times lower than expected. 
However a black body and adiabatic assumption is far from valid in such a hot star 

E2 NOTE TO EDITOR: compspec.ps should be over two columns 
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when one is considering amplitudes in visual. In a pure adiabatic assumption there is a 
180° phase difference between temperature and stellar radius (with maximum brightness 
coinciding with minimum radius). For the mode at 2.742 mHz (which do not appear to be 
contaminated by closely spaced modes), we find a phase difference of (-75±10)° implying 
that maximum brightness is slightly after maximum radius. More data are needed to clarify 
this result. 

EDITOR: PLACE FIGURE 3 HERE. 

E3 

There is evidence to suggest that the velocity amplitude is wavelength dependent. 
Figure 3 shows that the amplitude of oscillations at two frequencies decreases blueward of 
H(3. The amplitude of H8 is almost 25% smaller than that of H/3. If confirmed, this would 
be puzzling, given that all the Balmer lines are formed in the same part of the atmosphere. 

5. Conclusions 

The results presented here show that times-series spectroscopy of sdBs can be done 
using a medium-size (1.5 m) telescope. We have detected oscillations in velocity at 
previously published photometric frequencies. A longer times-series would allow most of the 
closely spaced modes in PG1605+072 to be resolved and thus help to confirm the variation 
of velocity amplitude with wavelength. Further studies may also determine whether the 
amplitudes of modes in this star are variable over time, or whether the apparent shift in 
primary pulsation mode is due to beating. We will also examine the behaviour of equivalent 
width and other line profile variations. 

E3 NOTE TO EDITOR: wave-vs-amp.dat.ps should be over one column 
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Table 1: Spectroscopic observations of PG1605+072. LS = La Silla; MS = Mt Stromlo. 



UT-date 


Site 


No. of 


No. of 






hours 


spectra 


1999 Julv 23-4 


LS 


4.01 


210 


1999 Tulv 24-5 


LS 


1.91 


105 


1 999 Tulv 25-6 


LS 


5 46 


280 


1999 July 26-7 


LS 


5.12 


280 


1999 July 30-1 


LS 


5.09 


260 


1999 July 31 


MS 


3.61 


168 


1999 July 31-1 


LS 


5.23 


280 


1999 Aug 01 


MS 


3.30 


155 


1999 Aug 01-2 


LS 


4.92 


280 


total 




38.65 


2018 
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Table 2: Photometric observations of PG1605+072, taken at the Sutherland site of the SAAO. 



UT-date 


No. of 


No. of 




hours 


obs. 


1999 Julv 20 


3.93 


666 


1 999 Tulv 21 


1 07 


1 82 


1999 July 23 


4.17 


679 


1999 July 24 


1.27 


211 


1999 July 28 


3.77 


597 


1999 July 29 


1.19 


199 


1999 July 31 


1.24 


210 


1999 Aug 01 


2.80 


454 


total 


19.44 


3198 
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Fig. 1. — Combined raw H/3 velocity and Johnson B light curve for PG1605+072. Velocity- 
points (in kms -1 ) are large dots while photometry points are small dots. Photometry in this 
Letter is presented in 'millimodulation units' or mmi, where the light curve is normalized by 
its own mean level (see Kilkenny et al. 1999). Observations were made on JD 2451380-393 
(see Tables 1 and 2). Note that positive velocity is redshifted by definition. 
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Fig. 2.— Amplitude spectrum of PG1605+072 in the 1000-5000 /iHz range. Panels on 
the left from top to bottom: R(3 velocity; H7 velocity (with no Stromlo data); B band 
photometry; Photometry by Kilkenny et al. (1999) in May, 1997 - the dotted lines represent 
the 10 strongest modes found. Panels on the right show greater detail centred around 
2100 /iHz. 
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Fig. 3. — Velocity amplitude vs wavelength for the 5 Balmer lines studied. The group of 
unresolved modes at 2.10 mHz is represented by diamonds, while the mode at 2.74 mHz is 
represented by triangles. 



